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Abstract Using data from the Heliospheric Imagers (His) onboard STEREO, 
it is possible to derive the direction of propagation of coronal mass ejections 
(CMEs) in addition to their speed with a variety of methods. For CMEs ob- 
served by both STEREO spacecraft, it is possible to derive their direction using 
simultaneous observations from the twin spacecraft and also, using observations 
from only one spacecraft with fitting methods. This makes it possible to test 
and compare different analyses techniques. In this article, we propose a new 
fitting method based on observations from one spacecraft, which we compare to 
the commonly used fitting method of Sheeley et al, (1999). We also compare 
the results from these two fitting methods with those from two stereoscopic 
methods, focusing on 12 CMEs observed simultaneously by the two STEREO 
spacecraft in 2008 and 2009. We find evidence that the fitting method of Sheeley 
et at, (1999) can result in significant errors in the determination of the CME 
direction when the CME propagates outside of 60° ±20° from the Sun-spacecraft 
line. We expect our new fitting method to be better adapted to the analysis of 
halo or limb CMEs with respect to the observing spacecraft. We also find some 
evidence that direct triangulation in the HI fields-of-view should only be applied 
to CMEs propagating approximatively towards Earth (±20° from the Sun-Earth 
line). Last, we address one of the possible sources of errors of fitting methods: 
the assumption of radial propagation. Using stereoscopic methods, we find that 
at least seven of the 12 studied CMEs had an heliospheric deflection of less than 
20° as they propagated in the HI ficlds-of-view, which, we believe, validates this 
approximation. 
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1. Introduction 

Thousands of coronal mass ejections (CMEs) have been observed remotely by 
coronagraphs and hundreds by in situ instruments since the 1970s (for a review 
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of CME observations, see, for example Hundhausen, 1993 and Howard, 2006). In 
the past five years, with the launch of spacecraft carrying heliospheric imagers 
(Coriolis and the Solar- Terrestrial Relations Observatory (STEREO)), we have 
witnessed the start of a new era, where CMEs can be routinely observed to 
radial distances as far as 0.5 AU with the Heliospheric Imagers (His, see Eyles 
et at, , 2009) and sometimes up to Earth's orbit. Their remote properties can 
now be compared to in situ measurements (Davis et al., 2009; Harrison et at, 
2009; Mostl et al., 2009). To analyze these new measurements, new methods 
and approximations must be devised. Shortly after the launch of Coriolis and 
STEREO, two simple ways to analyze wide angle heliospheric observations were 
developed: the Point-P (Vourlidas and Howard, 2006) and Fixed-$ (Kahler and 
Webb, 2007) approximations. Recently, another approximation has also been 
proposed (harmonic mean (HM), see Lugaz, Vourlidas, and Roussev, 2009) and 
these approximations have been completed by other techniques: visual fitting 
(Wood et at, 2009; Maloney, Gallagher, and McAteer, 2009), fitting to a family of 
pre-existing simulated ejections (Tappin and Howard, 2009; Howard and Tappin, 
2009) and fitting to known functions of the speed and direction (Rouillard et al., 
2008). These techniques can derive, in addition to the CME speed and position, 
its average direction of propagation. Two other methods have been proposed to 
analyze simultaneous (i.e. stereoscopic) CME observations in the heliosphere, 
by direct triangulation (Liu et al., 2010) or by a "tangent to a sphere" method 
(Lugaz et al., 2010), to obtain, at all times, the CME position and direction of 
propagation. 

Heliospheric instruments are planned in a number of future missions, including 
the Solar Orbiter and the Solar Probe. However, it is possible that there will not 
be other stereoscopic observations of CMEs by heliospheric imagers after the end 
of the STEREO mission. The time period from early 2008, when the STEREO 
spacecraft separation reached 45° to the early 2010 when it reached about 130° 
is the optimal period to have stereoscopic heliospheric observations. It should be 
used to validate, test and compare methods to analyze white-light heliospheric 
observations of CMEs. Such comparisons between methods was recently per- 
formed by Davis, Kennedy, and Davies (2010) between the visual fitting of COR 
images by Thernisien, Vourlidas, and Howard (2009) and the analytical fitting 
to a constant direction and velocity by Sheeley et al., (1999). The theoretical 
error and bias associated with the manual selection of the elongation data for the 
method of Sheeley et al., (1999) were also recently quantified in Williams et at, 
(2009). In the present article, we compare two different stereoscopic heliospheric 
methods with each other and with fitting methods to study the direction of 
propagation of CMEs. 

In this article, in Section 2, we propose a new fitting method, similar to 
that of Sheeley et al., (1999) but based on the model of Lugaz, Vourlidas, and 
Roussev (2009). We compare theoretically the two methods based on stereoscopic 
observations proposed by Liu et al., (2010) and Lugaz et al., (2010) in Section 3. 
In Section 4, we analyse in detail one of the 12 CMEs observed simultaneously 
by the two STEREO spacecraft in 2008 and 2009, before doing a statistical 
comparison of the methods based on these real analyses. The conclusions of this 
investigation are drawn in Section 5. 
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2. CME Direction of Propagation from Single- Spacecraft 
Observations 

In the ficld-of-view of an hcliosphcric imager, the simple assumption of "plane- 
of-the-sky" cannot be made to derive the position of a density structure. The 
position is therefore measured as the angle between the observing spacecraft, 
the Sun and the density structure, and it is commonly referred as the elongation 
angle. In the simplest case of a single plasma element or a spherical ejection, the 
elongation angle is a complex function of the heliocentric distance and of the 
angle of propagation with respect to the observing spacecraft. Therefore, the 
shape of the elongation vs. time curve depends on the speed and direction of 
propagation of a transient (Shceley et ai, 1999). When CMEs are observed to 
large elongation angles (up to 40° and beyond), the elongation vs. time profile 
can be fitted to analytical functions and the average speed and average direction 
of the CME can be derived under certain assumptions (see, for example, Webb 
et ai, , 2009 and Howard and Tappin, 2009). 

In the rest of this article, we note a as the elongation angle, and f3 as the 
direction of propagation of the CME, following the terminology of Rouillard 
et al., (2008). We also note dsT as the heliocentric position of the observing 
spacecraft (here STEREO) and t as the time. 



2.1. Methods 



We assume that an heliospheric imager observes a single plasma element, and we 
further assume that the direction of propagation is fixed. Under these assump- 
tions, the heliocentric distance, R, can be derived analytically as a function of 
the elongation angle: 

i?F$ = dsT- 



s'm(a + /3f$) 

This relation is usually referred to as the Fixed-$ (F$) approximation (e.g., 
Kahler and Webb, 2007). It can be inverted (Shecley et ai, 1999), assuming a 
constant velocity (CV), V = Rt, resulting in: 

/ Vp^cvt sin/3p$cv \ 

a = arctan . (1) 

\d S T - V F $cvt cos/? F $cW 

A measured profile of elongation vs. time can be fitted to a profile of calcu- 
lated elongations given by Equation (1). The Fixed-3> constant velocity (F$CV) 
procedure gives the average speed and direction of the transient. It was originally 
proposed by Sheeley et al, (1999) for slow ejections in the LASCO coronagraph 
fields-of-view. Rouillard et al, (2008) and Sheeley et al, (2008) applied this 
technique to corotating interaction regions observed by the STEREO/HIs. It 
has since been widely used to study CMEs (e.g., in Wood et al, , 2009, Davis 
et al, , 2009 and Davics et al., , 2009) in the hcliosphcric imager ficld-of-view. 

Lugaz, Vourlidas, and Roussev (2009) and Howard and Tappin (2009) have 
proposed a different way to derive heliocentric distances from elongation angles. 
In these works, the authors assume that the CME can be modeled as a circular 
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front whose center propagates on a fixed radial trajectory and which is anchored 
at the Sun. It is further assumed that the measured elongation angle simply 
corresponds to the angle between the Sun-spaceraft line and the line-of-sight 
tangent to this circular front. Here the direction, /3hm, corresponds to the di- 
rection of the nose of the CME, i.e. the direction of the point on the circular 
front with the largest heliocentric distance. As derived in Lugaz, Vourlidas, and 
Roussev (2009): 



This relation has been referred to as the harmonic mean (HM) approxima- 
tion, because the distance corresponds to the harmonic mean of the distance 
calculated with Fixed-$ and Point-P approximations. It can be also inverted, 
assuming a constant velocity (CV), as: 



A measured profile of elongation vs. time can also be fitted to profiles of calcu- 
lated elongations and the harmonic mean constant velocity (HMCV) procedure 
gives different estimates of the average speed and direction of the transient. 

2.2. Fitting Procedure 

In this section, we give a quick overview of the fitting procedure for the two 
methods just described. A more detailed explanation can be found in Rouillard 
et al., (2010) and Davis, Kennedy, and Davies (2010). First, a time-elongation 
map (J-map) is produced from a time sequence of running difference HI images. 
The procedure to produce J-maps is explained in details in Sheeley et al., (1999) 
and Davies et al., (2009) and we refer interested readers to these works. From 
a J-map, a number of sample points, N (typically 30 to 50), following the CME 
track are selected manually by an operator. These time series of time-elongation 
data are plotted on the Rutherford Appleton Laboratory website 1 and they were 
provided to us by C. Davis. The time-elongation data are the starting point of 
the current study. The typical error associated with the manual selection of the 
sample points has recently been addressed in Williams et al, (2009). They 
estimated the error in the direction to be typically 2-5° for CMEs observed up 
to 45° elongation and beyond. 

A time series is fitted with the theoretical formulae given by Equations (1) 
and (2) in the following manner. For a given value of the velocity, V, and the 



http: / /www. sstd.rl.ac.uk/stereo/HIEventList. html 



sm a 



1 + sin(a + /3hm) 




(2) 
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Elongation angle based on Fixed — Phi fit 



Elongation angle based on Harmonic — Mean fit 
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Figure 1. Time-elongation profile for a 400 km s — 1 transient propagating with different 
directions with respect to the observing spacecraft assuming the position is given by the 
Fixed-</> approximation (left) and the Harmonic Mean approximation (right). 

CME direction, /3, we calculate the standard error a (standard deviation of the 
residue) between the observed profile and the theoretical profile as follow: 



where ^observed is the observed elongation angle at time tk and ^theoretical is the 
elongation angle calculated with Equation (1) or (2) for the same time, tk is 
the sample time of the fc th point on the J-map track. We only use elongation 
measurements in the HI fields-of-view, because the speed of the CMEs in the 
coronagraphic ficlds-of-view is usually not yet constant. Therefore t\, the sample 
time corresponding to the first point on the J-map track is not obtained from 
measurements only. Solar measurements, such as flare time or first appearance 
in COR-1 field-of-view, can not be used to determine this time, since it typically 
corresponds to an elongation of 4 — 5°. For each value of V and j3, t\ is obtained 
by assuming a constant speed and by solving a t heoreticai(*i) = a b B erved(ii)- It 
yields a different value of t\ for each of the two fitting methods and for each 
value of V and (3. 

The procedure is repeated for values of the speed between 100 and 1000 
km s -1 by 1 km s^ 1 increment and for values of the direction between —30° and 
120° with 1° increment. This way, we obtain an "error map" giving the value 
of a for all possible combinations of V and (3. The best-fit values of (V, /?) is 
that for which a is minimum. We give the uncertainties in the fitting quantities 
corresponding to the value of (V, /?) for which a = 2<r m i n , corresponding to a 
95% certainty. A theoretical example is given in the following section, and an 
example based on actual data in Section 4. 

2.3. Theoretical Comparison 

In this section, we compare theoretically the two fitting methods. Figure 1 shows 
elongation vs time profiles given by Equations (1) and (2) for different directions 




(3) 



k=l 
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Figure 2. Left: Error in the fit for an analytical profile given by Equation (2) (HMCV) with 
a speed of 400 km s^ 1 and a direction of 40°. Color contours are the fitting error for the fit 
with Equation (1) (best fit marked by the black circle, blue is small error, white large) and line 
contours are the fitting error for the fit with Equation (2) (best fit marked by the white circle). 
Right: Best-fit direction (black diamonds) and velocity (red triangles) for a profile given by 
Equation (2) (HMCV) for a 400 km s CME observed from different directions and analyzed 
with Equation (f) (F0CV). 



of propagation. All cases correspond to a CME with a speed of 400 km s _1 
observed by STEREO- A. A similar figure was shown in Webb et ai, (2009) 
but the authors never used this as a fitting procedure to derive the direction of 
propagation of CMEs. 

As can be seen in this Figure, for a given direction, the new fitting method 
(HMCV) predicts that a CME exhibits less apparent deceleration/acceleration 
than as predicted by the F$CV fitting method of Sheeley et at, (1999). The 
reason is that, as a wide CME propagates close to the solar limb, an heliospheric 
imager observes more and more of the flank of the CME. It results in a apparent 
deflection towards the spacecraft of the observed structure and in a reduction of 
the "deceleration" due to geometrical effects. Conversely, for CMEs propagating 
towards the observing spacecraft, more and more of the CME nose is observed, 
resulting in less apparent acceleration. 

Next, we consider how a time-elongation profile given by Equation (2) is fitted 
with profiles from Equation (1) using the procedure explained in the previous 
section. This gives an idea of the theoretical error associated with the fact that 
the fitting method of Sheeley et ai, (1999) assumes that the heliospheric imager 
tracks a single plasma clement and not a front. 

In the left panel of Figure 2, we show an example of an error map correspond- 
ing to a CME front propagating 40° away from the Sun-spacccraft line with a 
speed of 400 km s _1 . A synthetic time-elongation profile for these values of V 
and (3 is generated using Equation (2). Fitting this profile with the procedure of 
Sheeley et at, (1999) results in a best-fit speed Vf$cv = 394 ± 6 km s _1 and a 
direction /?f$cv = 52° ± 4°. The error in the velocity is negligible but the error 
in the direction of propagation is relatively large. It is about twice as much as the 
error associated with the manual selection of points (Williams et al, 2009). This 
result holds true for other directions, as shown in the right panel of Figure 2. In 
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fact, we found simple relations between the direction of propagation and speed 
based on the F<3?CV fitting of Sheeley et ai, (1999) compared to that based on 
the HMCV fit: 



/Wv(°) = 2.067/Wv - 67.3 (4) 
AV{%) = 0.0168a| $cv - 2.12a F $cv + 66.38, (5) 

with AV = (Vhmcv — Vf^cv) /^hmcv and with correlations of 0.99985 and 
0.9988, respectively. This is found independently of the CME speed. It shows 
that around 60°, both fitting methods give the same direction. The F$CV fitting 
gives large errors in the direction for synthetic profiles corresponding to a wide 
CME with a direction greater than 80° or less than 40° with respect to the 
Sun-spacecraft line. 



3. CME Direction of Propagation from Stereoscopic Observations 

CMEs propagating between the STEREO-A and B spacecraft or wide CMEs 
propagating close to one of the STEREO spacecraft can be imaged to large 
elongation angles by the His onboard both STEREO spacecraft. Simultaneous 
measurements can be used to derive the CME direction of propagation for every 
pair of observations. Direct triangulation can be done under the assumption that 
both STEREO spacecraft observe the exact same plasma element (Liu et at, 
2010). Here, we derive a slightly different version of their formula, taking into 
account the difference in spacecraft heliocentric distances. The ratio of these 
distances is typically cLb/cLa = 1-07 and varies between about 1.04 and 1.14. 
Using the correct ratio results in a shift of the CME direction by as much as 
30° at large distances (beyond 0.5 AU) compared to assuming the ratio equal to 
one. 

The direction of propagation is given by 

/3Triang = arctan - ■ — - ■ '- , 6) 

\Pcos(a A + 7a) + cos(a B + j B ) J 



P = 



with 

ds sin as 
d,A sin a a 



where 7^4 and 75 are the separation between STEREO-A and B and Earth, 
respectively (both are defined as positive numbers). 

Alternatively, one can use the model of Lugaz et ai, (2010), which considers 
that the two STEREO spacecraft observe the tangent to a circular CME front 
anchored at the Sun. In this case, as shown in Lugaz et ai, (2010), the direction 
of propagation is given by: 



^Tang = /^Triang + aXCSm I — — — I , (7) 
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with 

Q = \/P 2 + 2Pcos{~f A +lB + a A + a B ) + 1, 

and /^Triang given by Equation (6). This Equation is correct for P cos(^a + &a) + 
cos(7b + Ob) > (which is equivalent to a a + ctB + 1A + 7b < n for = ds)- 
When the inequality is no longer true, the correct solution is: 

Arang = /^Triang - arcsin 

We refer to this method as the "tangent-to-a-sphere" method, because the elon- 
gation angle is understood as the direction of the tangent to the spherical CME 
front. 



4. Method Comparison Based on Real Data 

In this section, we compare these four methods on CMEs observed simultane- 
ously by the two STEREO spacecraft in 2008 and 2009. This allows us to use 
the stereoscopic methods and also to compare the fitting results for each of the 
two spacecraft. We found about 15 CMEs with stereoscopic observations in these 
two years, 12 of which had good enough data available on the RAL website for 
the study. These are the 2008 April 26, June 2, July 7, August 30 and December 
12 (two tracks: front and back) CMEs and the 2009 January 9, January 22, May 
9, May 13, September 4, October 18 and November 21 CMEs. Here, we analyze 
in details the 9 January 2009 CME, before comparing statistically the different 
methods in the next section. 

For an elongation vs. time profile, we judge one fitting method to be better 
than the other, when the best-fit error (as defined in Section 2.2) of this method 
is smaller than the best-fit error of the other method. It is also possible to 
compare the best-fit (V, (3) obtained with the data from STEREO-A with that 
obtained from the data from STEREO-B for the same CME. When there are 
values of (V, (3) within l-cr m i n of the best-fit value for STEREO-A and B data, 
we consider that there is good agreement between the two sets of data. A specific 
example is given in the following section. 

4.1. 2009 January 9 CME 

In January 2009, the STEREO spacecraft were almost in quadrature (separation 
of 89°). There was a weak CME on January 8 2009 around 20:00 UT detected 
in LASCO, COR2-A and COR2-B. It was observed in the His for about 100 
hours starting on January 9. A magnetic cloud was detected by STEREO-B 
starting on January 13 at 05:00 UT and ending at 22:00 UT the same day, and 
no magnetic cloud was observed by STEREO-A or ACE. From the predicted 
arrival time at 1 AU, we believe the January 9 2009 CME is a good candidate 
for this magnetic cloud. The January 9 CME was tracked until about 50° by 
STEREO-A (39 datapoints) and until about 60° by STEREO-B (52 datapoints). 
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Figure 3. Analysis of the 9 January 2009 CME. Top: STEREO-A data (cross) and best-fit 
solution obtained with HMCV (solid line) and F$CV (dash-dot line) fittings. Middle: Same but 
for STEREO-B data. Bottom: Combining the results of the HMCV fitting for both STEREO 
spacecraft. Line contours and color contours show the error for STEREO-A and B data, 
respectively (each contour is separated by rr m i n ). The pink striped area corresponds to value 
of V and f) simultaneously within the 95% certainty region of both fits. 

Wc fit the data from the two spacecraft with the two fitting methods. The 
time-elongation profiles corresponding to the best-fit values of (V, /3) are shown 
in Figure 3 for the two methods. The fitting with Equation (1) (F$CV) of the 
STEREO-A data yields a best-fit speed of 352 ± 61 km s -1 and a direction of 
-33.3° ± 17°. The fitting with Equation (2) (HMCV) of the same data yields a 
best-fit speed of 380 ± 75 km s" 1 and a direction of -55.3° ± 40°. The HMCV 
minimal error is about 22% larger than that the F$CV minimal error. The 
uncertainty for both methods is relatively large because the time-elongation 
profile is close to linear. As expected from our theoretical analysis, for this CME 
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Figure 4. Direction of propagation with respect to the Sun-Earth line of the 13 CME tracks 
from Table 1 with the tangent-to-a-sphere method (left) and the triangulation method (right). 
Note the different axis scales. 



observed almost at the limb by STEREO-A, the HMCV fitting yields a larger 
value of the direction with larger errors than the F$CV fitting. 

The F<J?CV fitting of the STEREO-B data yields a best-fit speed of 321 ± 10 
km s" 1 with a direction of -14.4° ± 7.5° for B data. The HMCV fitting of the 
same data yields a best-fit speed of 339 ± 6 km s _1 with a direction of —28.4° ± 
6° and the error is about 126% smaller than that from the F<3>CV best-fit. 
Because the time-elongation profile shows a strong acceleration starting around 
25°, the direction of the CME is much better constrained by the STEREO-B 
data than with the A data. STEREO-B observes the CME as a halo, and the 
result of a smaller direction based on the HMCV fit compared to the F$CV fit 
is consistent with our theoretical analysis. While the difference of 14° between 
the two directions may appear small, it is larger than the derived uncertainty. 
The much smaller error of the HMCV fit compared to the F$CV fit is a strong 
indication that STEREO-B did not observe the same part of the front at all 
time as assumed by the F$CV fitting. 

We can combine the results of the two spacecraft by looking for values of (V,/3) 
within l-cr m i n of the best values for STEREO-A and B data simultaneously. 
This is illustrated in the bottom panel of Figure 3 for the HMCV fit and it 
greatly reduces the uncertainties in the best-fit values. Combining the two set 
of observations, the best-fit value corresponds to a CME with a speed of 339 ± 
6 km s" 1 and a direction of —30° ± 4.5° (pink striped region). We perform the 
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Table 1. Direction of propagation for the 13 CME tracks from the four methods: F<I>CV fitting method of Sheeley et al.l 
(1999) for STEREO-A and STEREO- A data, HMCV fitting method for STEREO- A and STEREO-B data, triangulation 
and tangent to a sphere methods for columns 2 to 7, respectively. For the fitting methods, bold results indicate that this 
method has the lowest error result for this dataset. When there are in situ measurements, we indicate it in column 8 as 
well as the spacecraft angular separation with Earth, (b) refers to the back of the December 12 CME. 



CME 


F*CV A 


F<t>CV B 


HMCV A 


HMCV B 


Triang. 


Tangent 


in situ 


2008 
















Apr. 26 


-33.5° ±18 


2.1° ± 7 


- 34.5° ± 30 


-19.9° ± 10 


-9° ± 4 


-34° ± 5 


B (-24°) 


June 2 


-24.2°± 6 


20.9°± 12 


-34.2° ± 8 


-0.1° ± 25 


-3.5°± 1 


-16.5°± 3 


B (-25°) 


July 8 


-24.6°± 7 


2.6°± 11 


-34.6° ± 8 


-25.4° ± 21 


-3.5°± 2 


-17°± 6 


none 


Aug. 30 


-0.4°± 8 


19.2°± 11 


-19.4°± 25 


12.2°± 27 


4°± 5 


5°± 14 


none 


Dec. 12 


-14.7° ± 13 


12.6° ± 7 


-14. 7° ± 19 


18.6°± 11 


3°± 4 


8°± 10 


ACE 


Dec. 12 (b) 


8.3°± 5 


-1.5°± 7 


25. 3° ± 13 


-14.5°± 11 


8.5°± 3 


20° ± 6 


ACE 


2009 
















Jan. 9 


-33.3°± 17 


_14.4°± 8 


-55.3°± 40 


-28.4°± 6 


-6°± 4 


-11.5°± 10 


B (-46°) 


Jan. 22 


-10.1°± 19 


22°± 8 


-4.1°± 30 


34° ± 20 


7°± 6 


17° ± 12 


none 


May 9 


10.4°± 6 


8.1°± 12 


18. 4° ± 9 


5.1°± 32 


12°± 4 


23°± 11 


none 


May 13 


-14. 3° ± 8 


18.1°± 10 


-27.3° ± 16 


23.1°± 22 


12.5°± 8 


16°± 25 


none 


Sep. 4 


5.6°± 11 


-19. 5° ± 4 


15.6° ± 16 


-32.5°± 13 


-18°± 4 


-50°± 11 


none 


Oct. 18 


-19. 8° ± 12 


-16.1°± 7 


-43. 8° ± 26 


-23.1°± 11 


-18°± 8 


-30°± 17 


B (-60°) 


Nov. 21 


-18.7°± 10 


-22. 6° ± 6 


-44. 7° ± 17 


-29.6°± 14 


-15.5°± 4 


-30° ± 4 


unknown 



same procedure for the F<I>CV fit (not shown here) and we find a best-fit value of 
the speed of 319 ± 8 km s _1 with a direction of —18.5° ± 2.5°. If we assume that 
the magnetic cloud observed at STEREO-B on January 13 and not observed by 
ACE is associated with this ejection, and taking into account that STEREO-B 
was 46.6° away from Earth, the direction of -30° ± 4.5° from the HMCV fitting 
method appears more likely to be the correct one that —18.5° ± 2.5° from the 
F$CV fitting. 

Applying the stereoscopic methods to this CME, we find a direction of — 13° ± 
7° with the tangent-to-a-sphere method and —6.5° ± 3° for the triangulation 
technique limiting our analysis to distances between 30 and 120 Rq. The vari- 
ation with distance of the direction of propagation from these two stereoscopic 
methods is shown in Figure 4. Both stereoscopic methods give results not in good 
agreement with the fitting methods and with a significant deflection towards the 
East. One interpretation is that the fitting methods give the direction of propa- 
gation at large heliosphcric distances. If we assume a linear trend with distance, 
triangulation predicts a direction of about —15° at 1 AU and the tangent-to-a- 
spherc method predicts about —30°. Based on these values, there is relatively 
good agreement between the two methods based on the Fixed-;/) approximation 
(F$CV and triangulation) and also good agreement between the two methods 
based on the harmonic mean approximation (HMCV and tangent-to-a-sphere). 
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Direction RAL Best Spacecraft (°) Direction New Fit Best Spocecroft (°) 

Figure 5. Comparison of the CME direction of propagation with respect to the Sun-Earth 
line obtained from fitting the observations by STEREO-A with that obtained from fitting the 
observations by STEREO-B with Equation (1) (F*CV, left) and Equation (2) (HMCV, right). 

4.2. Summary of the Analysis 

The full results of our analysis are shown in Table 1, where directions are given 
with respect to the Sun-Earth line (positive number is towards the West, i.e. 
in the direction of STEREO-A). In this Table, the results of the fitting method 
with the smallest error are marked in bold (better fit). In addition to the 9 
January 2009 CME, we are able to match the results of both fitting methods for 
STEREO-A and B spacecraft for the 2009 October 18 and November 21 CMEs. 
The results for the back of the 12 December 2008 CME can be similarly matched 
for the F$CV fitting method of Sheeley et al, (1999), but not with the HMCV 
fitting method. The contrary is true the 8 July 2008 CME (can be combined 
with the HMCV method but not with the F$CV method). 

In addition to this analysis of CMEs observed by both STEREO spacecraft, we 
also analyze four CMEs which are observed away from the Thomson sphere: the 
23 May 2008 CME observed by STEREO-B with a direction of 13° with respect 
to the Sun-spacecraft line, the 16 August 2007 CME observed by STEREO-A 
with a direction of 32° away from the spacecraft, the December 4, 2007 CME 
observed by STEREO-B with a direction of 75° away from the spacecraft and 
the March 18 2008 CME observed by STEREO-A with a direction of 80° away 
from the spacecraft. The results of the analysis of these four CMEs observed by 
only one STEREO spacecraft are shown in Table 2. 

4.3. Comparing the Fitting Methods 

As shown in Table 1, the HMCV fitting method results in a better fit (as defined 
by a smaller fitting error) to the data about half the time (12/26 with one in- 
stance where there is no improvement). Comparing the direction of propagation 
obtained from STEREO-A data to the direction obtained from STEREO-B data, 
the HMCV fitting method yields results in slightly better agreement (see Fig- 
ure 5), although statistically the cross-correlation between the direction obtained 
from STEREO-A and the direction obtained from STEREO-B only increases to 
0.214 from —0.125. To detect any bias in the methods, it is enlightening to look 
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at the direction of propagation with respect to the observing spacecraft. For 
this analysis, we treat all the observations as individual events and include the 
four CMEs from Table 2, resulting in 30 datapoints. The results are shown in 
Figure 6, where better fits with the new method are shown with black squares 
and worst fits with red triangles. From this Figure, it appears clear that there 
is a systematic bias associated with one of the fitting methods, although it is 
slightly less pronounced than what we found in the theoretical analysis (Section 
2.3). In fact, we can find a linear relation between the direction derived from the 
two methods with a cross-correlation of 0.966: 

/W°) = l-835/3 F$ -45. (8) 

This relation is valid when we consider only the better fits or the entire dataset. 
This shows that the two methods result in CME directions within ±15° of each 
other for direction of 55° ± 18°, similar to what we found in the theoretical 
analysis. It is worth noting that out of these 30 CME tracks, the direction given 
by the F<f>CV fitting and the direction given by the HMCV fitting are within 
each other 95% confidence range in all but ten cases. All these cases have a CME 
predicted direction of more than 70° or less than 20°. For only one of these ten 
cases, the error of the F<3?CV fitting is less than that of the HMCV fitting (plus 
one case where they are identical). 

Out of these 12 CMEs, at least five resulted in a magnetic cloud observed in 
situ by one of the spacecraft, providing some information of the real direction 
of propagation. Considering only the best observing spacecraft (the one which 
observes the CME to the largest elongation angle), both methods would have 
successfully predicted a CME hit with two exceptions: the 12 December 2008 
CME, for which the HMCV fitting method may have predicted a miss at ACE or 
a hit also at STEREO-B and the 18 October 2009, for which the F$CV fitting 
method would have predicted a miss at STEREO-B. 

Five of the studied CMEs were also part of the study by Thernisien, Vourlidas, 
and Howard (2009). They are the 4 December 2007 CME with a direction of 71°, 
2008 March 18, April 26, June 2 and July 8 with direction of -83°, -21°, -37° 
and —23°, respectively. For three of them, the direction obtained from the HMCV 
fitting method is closer to this direction from COR data than is the direction 
from the F$CV fitting method. For the two other (April 26 and July 8), it is 
comparable or better if one uses the best-fit value from the spacecraft observing 



Table 2. Direction of propagation for the four additional CMEs, selected because of the small or largi 
direction of propagation with respect to the observing spacecraft. 





CME 


F$CV 


HMCV 


Remarks 


2007 Aug. 16 

2007 Dec. 4 

2008 Mar. 18 

2008 May 23 


-17.6° ±9 
53.6°± 7 
-56.9°± 12 
-11.7°± 9 


8.4° ± 18 
69.6°± 10 
-84.9°± 16 
-37.7°± 4 


Same fitting error 
71° from Thernisien, Vourlidas, and Howard (2009) 
—83° from Thernisien, Vourlidas, and Howard (2009) 
Fit improved by 160% 
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Direction RAL (°) Direction Tangent 

Figure 6. Left: Comparison of the best-fit direction of propagation with respect to the Sun- 
spacecraft line obtained from the F<I>CV fitting with that from the HMCV fitting. Red and 
black data points represent values for which the best-fit is obtained with the F<3?CV and 
HMCV fittings, respectively and the green points are for cases when the best-fit error with 
the two methods is the same. Right: Comparison of the average direction obtained from the 
triangulation and from tangent-to-a-sphere method. 

the CME as a halo. In fact, it is clear from our study that the direction is best 
constrained when it is less than 35° or more than 75°, but when the CME can be 
tracked over a long-enough elongation range. Taking into account our theoretical 
analysis, the comparison of the direction with in situ measurements and with 
the procedure of Thcrnisicn, Vourlidas, and Howard (2009), we conclude that 
the fitting method of Sheeley et ai, (1999) has a systematic bias for CMEs 
propagating beyond 60° ± 20° away from the Sun-spacecraft line. 

4.4. Comparing the Stereoscopic Methods and CME Heliospheric Deflection 

Figure 6 shows the comparison of the two stereoscopic methods with each other. 
It also appears that there is a systematic bias, the triangulation method only 
yields directions of propagation very close to the Sun-Earth line. In fact, the two 
terms in Equation (7) are almost always of the same order so that /?Tang ~ 2 x 
/?Triang- It is less clear which of the two methods is systematically biased. Based 
on observations at 1 AU, triangulation works best for these CMEs propagating 
close to the Sun-Earth line such as the 12 December 2008 CME and the 4 
September 2009 CME. It is also, in general, much less noisy than the tangent- 
to-a-sphere method. However, it gives inconsistent results for CMEs observed 
as a halo by one of the spacecraft (26 April and 2 June 2008 and 9 January 
2009). For these CMEs, the stereoscopic method of Lugaz et ai, (2010) or a 
more complex analysis (Wood and Howard, 2009; Tappin and Howard, 2009) 
must be used. This result is expected as triangulation requires both spacecraft 
to observe the same plasma element. However, it points to a strong limitation of 
the triangulation method: it can only be used if the direction can be estimated to 
be close to 0° by another method. To use triangulation in the HI ficld-of- views, 
it is necessary but not sufficient that both STEREO spacecraft observe the CME 
to large elongation angles. 

Next, we study the assumption of constant heliospheric direction. We focus on 
the tangent-to-a-sphcre method, since the derived directions and their variations 
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are about a factor of two lower based on the direct triangulation. This provides a 
worst-case scenario. For all of the 2008 CMEs except the August 30 ejection (four 
CMEs), the results are consistent with a very small heliospheric deflection (less 
than 15° ALT 1 in absolute value). This also true for the 21 November 2009 CME, 
if we only consider the measurements in HI-2 field of view (from 0.5 to 1 AU). 
For 2 other CMEs (2009 May 9 and September 4), the heliospheric deflection is 
found to be about 30° AU -1 , while the other five CMEs are found to have larger 
heliospheric deflection or noisy data. This difference between CMEs observed in 
2008 and those observed in 2009 may indicate that stereoscopic measurements 
can only be made for moderate spacecraft separations (less than 90°). Taking 
into account the measurement errors and the limitations of the stereoscopic 
methods, we believe that a deflection by less than 20° during the duration of the 
observations, is consistent with the assumption of radial propagation in the HI 
fields-of-view. This is the case for seven of the 12 studied CMEs. 

5. Discussions and Conclusions 

In this article, we have examined some of the methods used to derive the 
direction of propagation of coronal mass ejections (CMEs) from observations 
by heliospheric imagers. We have focused on four methods: the fitting method 
of Sheclcy et al, (1999), the triangulation method of Liu et al, (2010), the 
tangent-to-a-sphere method of Lugaz et al, (2010) and a new fitting method 
derived in this article. 

Because it is impossible to know the exact direction of propagation of a CME, 
we have focused our analysis of real data on ejections observed simultaneously 
by the two STEREO spacecraft. For these CMEs, it is possible to use the 
stereoscopic methods and also to compare the direction of propagation based 
on fitting methods for each of the instruments. Overall, we have found a very 
small or no correlation between the direction obtained from observations by one 
spacecraft and that obtained from observations by the other spacecraft, and this, 
both for the fitting method of Sheeley et al, (1999) and for the new method we 
proposed. However, for four out of 13 CME tracks, it is possible to combine the 
results of the fit based on the data from each spacecraft within the 95% certainty 
interval to find a well-defined direction. 

We have addressed two of the three possible sources of error for the fitting 
method of Sheeley et al, (1999) applied to CMEs: i) the assumption of radial 
propagation (absence of heliospheric deflection), and ii) the assumption that 
the CME is of negligible width (or that the same part of the CME is always 
observed). We addressed i) by using stereoscopic methods which can give at all 
time the CME direction of propagation instead of its average value. Stereoscopic 
methods need to be better tested and validated but for about nine out 13 CMEs, 
we have found reasonable results with the method of Lugaz et al, (2010) and 
for two CMEs with the triangulation method of Liu et al, (2010). For the other 
two cases, the method of Lugaz et al, (2010) is too noisy and triangulation 
results in a large underestimation of the CME direction of propagation. In seven 
of these good cases, we have found no sign of strong heliospheric deflection of 
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CMEs, the total deflection being less than 20° in the HI field-of-view. This result 
validates the usage of fitting methods to determine CME direction, since these 
methods assume no heliospheric deflection. 

We have proposed a new fitting method which takes into account the CME 
width. We have found that the assumption of negligible CME width yields worst 
fit for the method of Sheeley et ai, (1999) compared to this new method for 
CMEs propagating outside of 60° ± 20° from the Sun-spacecraft line. We believe 
this new method should be used for halo and limb CMEs as seen by the observing 
spacecraft. However, both methods can be easily used for all observations and 
the one with the smallest residual error should be chosen to provide the CME 
estimated speed and direction on a case-by-case basis. In our analysis of 30 CME 
tracks, we have found that each fitting method gives a better fit in about half 
the cases. We have also found that, when a CME is observed to large elongation 
angles by both STEREO spacecraft, the CME direction is best determined using 
the data from the spacecraft which observed the CME as a halo. This is especially 
true with the new fitting method. 

There is a third issue with the two fitting methods: the assumption of constant 
heliospheric velocity. All the CMEs in our sample (with the exception of the 26 
April 2008 with a speed of about 600-700 km s _1 ) have speed between 280 
and 450 km s _1 and we believe this assumption is approximatively true. This 
issue is hard to address because of the lack of fast CMEs in the past three 
years. Numerical simulations combined with a proper treatment of the Thomson 
scattering (Lugaz et ai, 2008; Lugaz et ai, 2009; Manchester et ai, 2008; Odstrcil 
and Pizzo, 2009) are one of the only ways to address this question of the effect 
of CME deceleration and we plan to look into this in future research. 
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